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ABSTRACT 

The variability of the continuum spectral energy distribution has been 
analyzed for a complete magnitude-limited sample of quasars in Selected Area 
57, observed at two epochs in the photographic U, Bj, F and iV bands with the 
Mayall 4m telescope at Kitt Peak. Changes 5a of the spectral slope a appear 
correlated with brightness variations 5logf u indicating an average hardening 
of the spectrum in the bright phases. This confirms that the correlation of 
variability with redshift, found in a single observing band, is due to intrinsic 
spectral changes. The average observed 5a — 5logf v relation is consistent with 
the spectral change due to temperature variation of a black body of about 
2.5 • l^K. 

Subject headings: galaxies: active - galaxies: multicolor photometry - galaxies: 
Seyfert - photometry : variability - quasars: general 
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1. Introduction 

Variability provides a powerful tool for constraining the physics of Active Galactic 
Nuclei. The original models based on a central black hole surrounded by an accretion disk 
( [Rccs 19841) were based on the constraints on the source size and energy density provided 
by the early single-band variability studies. Single-epoch multi-wavelength observations, 
covering simultaneously the range from radio to gamma-ray frequencies, are of crucial 
importance in suggesting the physical mechanisms responsible for the emission. However, 
in general they cannot sufficiently constrain the large number of parameters involved in 
models. So far mult i- wavelength monitoring, with adequate time sampling and duration, 
has been possible for a small number of low-redshift objects, thanks to large international 
collaborations (see the review of |Ulrich, Maraschi, and Urry 19971 an d refs. therein). These 



observations provided further clues to and constraints on the radiation processes and 
reprocessing mechanisms, since changes of the spectral energy distribution (SED) associated 
with flux variations can be interpreted in terms of an interplay of emission components of 
different spectral shape and variability properties. Various classes of active galactic nuclei 
(AGN) show different variability properties. Blazars show strong (> 1 mag) variations of 
the observed fluxes on time scales from days to months from radio to gamma frequencies, 
while in the optical-UV Seyfert 1 and normal quasars (QSOs) show smaller (< 0.5 mag) 
variability on time scales greater than a few months, although > 1 mag variations on time 
scales of days have been detected in X-rays in some Seyfert galaxies (Mushotzki et al. 
1993 and refs. therein) (in the following we ignore the traditional distinction at M#=-23 
between Seyfert 1 galaxies and QSOs). The result of multi- frequency monitoring indicates 
that the optical-UV continuum of AGNs hardens during the bright phases flCutri et al. 



1985| , Edelson et al. 1990| , [Kinney et al. 1991 , |Paltani and Courvoisier 1994| ) in the near 



infrared-optical-UV range, while little or no correlation is thought to be present in BL Lac 
objects flEdelson 1992| , see however [Massaro et al. 1996| ). 
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Optical observations at faint magnitudes (B pa 20 — 23 mag) enable magnitude-limited 
statistical samples to be constructed (mostly of radio-quiet objects) in a single field of a 
large telescope, and repeated broad-band imaging provides light curves for each AGN in 
the sample ([Hawkins 1983| , |Koo fc Kron 19881 , |Trevese et al. 1989| (T89), [Trevese et al. 



|19941 (T94), ICristiani, Vio, fc Andreani 1990| , |Hawkins fc Veron 1993| , |Hawkins fc Veron 



|1996| , |Bershady, Trevese fc Kron 1998| (BTK98)) allowing the ensemble variability properties 
of the sample to be assessed. 

Studies of this kind are complicated by the fact that in all magnitude-limited samples 
the luminosity is strongly correlated with redshift (L-z correlation), since the majority 
of objects lie close to the limiting flux. This makes it difficult to isolate the intrinsic 
variability-luminosity (v-L) and variability-redshift (v-z) correlations. 

Since early studies, the relation between absolute optical luminosity and optical 
variability has been addressed by several authors ( [Angione fc Smith 1972| , |Uomoto et al. 



|1976| , |Pica fc Smith 1983| ), with the aim of discovering whether QSOs are made up of 



multiple sources or single coherent sources. In the former decrease of variability 

amplitude with luminosity is expected. In recent years, consensus has progressively grown 
that brighter QSOs are less variable on average flTrevese et al. 1994| (T94) JHook et al. 



1991 ICristiani et al. 19961) (see Giveon et al. 19991 for a recent summary). This would seem 



to support a model where the active nucleus is powered by a series of supernova explosions 
Aretxaga fc Terlavich| , [Aretxaga, (Jid Fernandes fc Terlavich| , though this model can fail 



to reproduce the variability amplitude in the case of most luminous QSOs, according 
to [Hawkins 2000| . A similar negative correlation of X-ray variability with the 1-10 keV 
luminosity was found in Seyfert galaxies and QSOs and interpreted in terms of large 
number of incoherent flaring subunits by |Green et al. 1993| . On the contrary, a positive v-L 
correlation was found by [Edelson 1992j for the Blazar population. 
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Concerning the L-z correlation, the situation is further complicated by the fact that 
the different variability indicators, used by different authors, are affected in different ways 
by the v-L and L-z correlation and by spurious redshift dependence due to the combination 
of time dilation, intrinsic variability time-scales, and total duration of the observational 
campaign. This led Giallongo, Trevese & Vagnetti (1991) (GTV) to introduce a variability 
indicator based on the rest-frame structure function. Adopting this indicator, they found 
a positive correlation of variability with redshift, taking properly into account the effect of 
v-L and L-z correlations. This result has been subsequently confirmed by the analysis of 
the QSO structure function in bins of luminosity and redshift, performed by |Cristiani et al 
1996| using all the statistical samples available at that time. 



GTV suggested that the v-z correlation can by explained by a hardening of the 
spectrum in the bright phases (and vice versa), coupled with the increase of the rest-frame 
frequency with redshift, for a fixed observing band. 



Di Giemente et al. 1996, on the basis of the analysis of various QSO samples, have 
shown that indeed, on average, variability increases with the rest-frame frequency. This 
provides an indirect statistical argument for a slope change of the quasar SED associated 
with luminosity variability. Moreover the increase of variability with rest-frame frequency 
is quantitatively consistent with the v-z correlation. 

However, a direct statistical quantification of the change of spectral slope associated 
with luminosity variation is still missing. Moreover, even though variability has historically 
played a key role in constraining models of the QSO central engine, the physics of variability 
is still largely unknown. In fact, quite diverse variability mechanisms have been proposed, 
including supernovae explosions ( [Aretxaga fc Terlavich| ) , instabilities in the accretion disk 



( [Kawaguchi et al. 1998| ), and gravitational lensing due to intervening matter ( [Hawkins fc 



|Veron 1993| , |Hawkins 2000| ). Discriminating between these different models on the sole basis 
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of single band correlation functions and v-L and v-z correlations is a very difficult task. 

In the present paper we analyze UBjFN photometry of the faint QSO sample of SA57 
to look for direct evidence of an average spectral hardening for increasing flux. A statistical 
quantification of the spectral slope changes as a function of luminosity variations provides 
a new constraint for models of the variability mechanism. We also discuss the consistency 
of the SED variations with the v-z correlation and with temperature changes of an emitting 
black body. 

This paper is organized as follows: data and reduction procedures are described in 
section 2; the statistical analysis of SED changes and their relation with flux variations are 
described in section 3; a concluding discussion is presented in section 4. 



2. The AGN Sample and the Spectral Energy Distributions 

Statistical studies of variability require a QSO sample based on clearly known selection 
criteria, whose effects are clearly understood and quantifiable. Among the existing QSO 
samples, the magnitude-limited sample of faint QSOs of Selected Area 57 ( [Kron fc Chiu| 



|1981| , [Koo, Kron fc Cudworth 1986| (KKC)) is probably the most studied. This field has 
been repeatedly observed since 1974 with the Mayall 4 m telescope at Kitt Peak in the 
photographic U, Bj, F, N bands. Different techniques, such as colors, lack of proper motion, 
and variability, have been applied to search for AGNs down to Bj = 23 and to estimate the 
completeness of the sample. Spectroscopic confirmation of almost all of the QSO candidates 
has been obtained (|Kron fc Chiu 1981| , |Koo, Kron fc Cudworth 1986| (KKC), |Koo fc Kron 



|1988| ,T89, |Majewski et al. 1991| ,T94), with the exception of the faint part of the sample of 
extended sources of BTK98. This allows for the first time to perform a statistical analysis 
of the variability of the spectral energy distribution of QSOs. Single-band variability 
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studies (T89, T94, BTK98) of the AGN sample where based on the entire collection of J 
(or Bj) plates. However, as discussed in the next section, the most reliable analysis of 
SED variability is obtained considering only "simultaneous" U, Bj, F, N observations. For 
this reason we selected the two sets of plates taken in April 1984 and in April 1985. The 
list of the plates is given in Table 1, which contains plate identification, the date, and the 
photometric band. The effective bandwidths of the U, Bj, F, N bands are 570, 1250, 1150, 
1720 A respectively and the effective wavelengths are 3540, 4460, 5980, 7830 Arespectively. 

EDITOR: PLACE TABLE g HERE. 

Photometric methods and signal-to- noise ratio optimization are described in T89, T94 
and BTK98. The AGN sample consists of 40 spectroscopically confirmed objects. Table 2 
gives the photometry on each of the 9 plates considered in the present analysis. 

EDITOR: PLACE TABLE ?? HERE. 

Of these 40 objects, 35 are the QSOs discussed in T94, while another 5 objects are 
the spectroscopically confirmed AGNs selected by BTK as variable objects with extended 
images. The limiting magnitude of the sample is ni2 < 23.5 mag, where rri2 is the Bj 
magnitude evaluated in a fixed aperture of radius 0".5 (but see BTK98 for faint extended 
objects). The cumulative redshift distribution has 25%, 50%, and 75% of the objects below 
z=(0.75, 1.16, 1.80) respectively. 

EDITOR: PLACE FIGURE HERE. 

Figure 1 shows the Hubble diagram of the sample. Figure 2 shows a comparison of the 
U,Bj,F fluxes at one epoch (April 1984), normalized at A = 2000 A, as a function of the 
rest-frame frequency. 
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EDITOR: PLACE FIGURE | HERE. 

Most objects have similar spectral slopes, with five main exceptions, which are objects 
with a much steeper spectrum. Two of them are low redshift and two are the highest 
redshift objects. The former two, N ser 110195 with z=0.243 and N ser 114264 with z=0.287, 
belong to the sample of BTK98 and are relatively faint AGNs, with M# 7 = —19.6 and 
Mbj = —20.7 respectively. Their r.m.s. variability is among the lowest in the sample of 
BTK. Since we know that, in general, fainter AGNs show a stronger variability, we may 
argue that in this case the light from the nucleus is diluted by the stellar light of the host 
galaxy. This could also explain the steep spectrum. If so, the spectral changes would also be 
affected by the presence of a constant stellar component. Two other steep-spectrum objects 
are the most distant in the sample, N ser 111610 with z=3.02 and N ser 101392 with z=3.08, 
and Mbj = —20.7 and Mb t = —24.9 respectively (bottom right in Figure 2). Their spectral 
slopes are strongly affected by intergalactic Ly-a absorption. In particular the Bj band 
falls in the (rest-frame) region between the emission Ly-a and the Lyman-limit, and the U 
band falls beyond the Lyman-limit. The relevant fractional attenuation of the continuum 
can be evaluated as ^ 0.2 and and ^ 0.5 respectively, e.g. from |Steidel fc Sargent (1987) 



Notice that, in Figure 2, the downwards shift of the F-band flux is due to the normalization 
at 2000 A computed from the extrapolation of the steep uncorrected spectrum, and is 
dramatically reduced if the extrapolation to 2000 A is done after the correction for the 
Ly-a absorption. Once corrected for the absorption the two spectra become consistent 
with the average slope (a) ~ —0.5 of the sample (see |Gaillongo, Gratton fc Trevese 1990|) . 
We stress that the changes of the spectral index due to variability are independent of the 
intergalactic absorption, even though the apparent SED is affected. Figure 3 shows, with 
arbitrary scale, the data reported in Table 2. 



EDITOR: PLACE FIGURE | HERE. 
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In the majority of the cases the deviations of the SED from a linear logf u — logv 
relation are small, but some objects have a local maximum of emission in the observed 
frequency range. Of course, such a humped SED can be relevant for physical models of 
the emission mechanism. Correlations of colors with other properties, e.g. the X-optical 
slope (|Miyaji et al. 19971 ), can be affected by artificial curvature due to non- simultaneous 
measures of different bands. 



3. Variability of the continuum spectral energy distribution 

Our aim is to correlate the SED changes with the brightness variations. We do this 
by correlating the variations 5a of the spectral slope a, defined by logf u oc u a , with the 
variation of 5logf u , where v is the rest-frame frequency of one of the observing bands. In 
order to minimize the effect of noise, we are subject to several restrictions. On average the 
best signal-to-noise ratio is obtained on J plates; thus we use pairs of J plates to compute 
Slogfvj. The noise in a is lower for larger leverage on the logv axis and the larger the 
number of bands used. The data used for computing 5logf Uj must not be used in the 
calculation of a, otherwise a spurious correlation is found, namely 5a and 5logf Uj appear 
correlated even in the absence of variability, due to the presence of correlated noise in 
both 5a and 5logf Vj . This statistical bias may even lead to strong and highly significant 
correlation ( [Massaro fc Trevese 19961) . It must be considered also that the SED cannot 



be represented exactly by a power law in the wavelength range studied. This is not a 
problem as long as a at different epochs is computed in a similar way. In fact a is a linear 
combination of the values of logf Uj measured in different bands, thus any deviation of the 
SED from a power law, while affecting a, does not affect 5a, which is a linear combination 
of 5logf Vj values. This is true even when the deviation from the power law is produced 
by a (non-variable) absorption in some of the bands. However, computing a in a very 
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different way at two epochs, e.g. using U and iV or J and F bands, can generate 8a values 
not due to variability. This 'spurious' 8a can even be different for two identical QSOs, if 
they have different redshifts. The above considerations led us to select the two-epoch set 
of plates described in Table 1, to minimize spurious effects and to increase the reliability 
and robustness of the results. We determined a from a linear fit in the logf v — logv plane 
from U, Bj, F, N data at the first epoch and from U, F, N data at the second epoch (the 
N magnitude of N ser 104855 is missing at the first epoch). In this way two more J plates 
(MPF3919 and MPF3977), one at the first and one at the second epoch, are available to 
compute an independent 8f Uj . Emission lines, which are known to vary on the same time 
scale of the continuum emission, may also affect the spectral slope variations. Whether 
they cause a positive or negative contribution to the slope depends on the particular band 
and redshift. Of course, if line variation were perfectly synchronous and proportional to 
continuum variation they would not cause any additional slope variability. A recent paper 
of [Kaspi et al. 2000], provides light curves of H a , Hp, H 1 lines and continuum of 28 PG 



QSOs with a sampling interval of 1-4 months for a total observing period of 7.5 years. 
From the line-continuum cross-correlation function the authors derive time delays of the 
order of 100 days for 17 of the 28 QSOs. H a flux fluctuations are smaller than 20% of the 
continuum fluctuations. Continuum and line variation are, in general, strongly correlated 
and the time width of the cross correlation function, which depends on the power spectrum 
of variability, is larger than the delay. As a result, in most cases the cross-correlation at 
zero time-lag is large 0.5) so that line variations synchronous and proportional to the 
continuum variation can be considered a good approximation for the purpose of the present 
analysis. 

Although we cannot correlate the slope changes with flux variations for each individual 
object, we can study the "ensemble average" of spectral variations of the 40 objects of our 
sample between two epochs. In Figure 4, 8a versus 8logf Vj is reported for the entire sample, 
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with the two linear regression lines, showing an average increase of the spectral slope a for 
increasing fj, i.e. a hardening of the spectrum in the brighter phases. The distribution of 
5a's is symmetric about zero, indicating the absence of any strong systematic effect on the 
spectral slopes. 

EDITOR: PLACE FIGURE | HERE. 

The correlation coefficient between 5a and Slogf Uj is p = 0.39 and the probability of the 
null hypothesis is P(> p) = 1.28 ■ 10~ 2 . We take the above result as statistical evidence for 
a relation of the type 5a = a + b 5logf Uj between a and luminosity variations. To check 
the dependence of the result on the two most deviant points (iV ser .=118122 Alogfj=-0.276, 
and iV se7 .=104885 Alogfj=0.292) we evaluated the correlation excluding in turn: i) both 
the points, ii) the first point and iii) the second point. We obtained respectively: i) 
p = 0.288 P{> p) = 7.98 • 10- 2 , ii) p = 0.290 P{> p) = 7.33 • 10~ 2 and iii) p = 0.395 
P(> p) = 1.29 ■ 10~ 2 . To determine the slope b of the above relation, the errors on 
5a and 5logf Vj must be taken into account, especially if they differ on the two axes. 
Errors on 5logf v for the plate pairs of each band can be computed as follows. From the 
photometric catalog of the area of sky which has been monitored (KKC,T89,T94,BTK98), 
each containing 2540 objects brighter than Bj=22.5, we computed 5logf Uj for each object 
and for each band. Then, in each band, we sorted the objects according to their magnitude, 
and for each object we took the 100 nearest neighbors (in magnitude) and computed the 
r.m.s. value of 5logf Vi (i=U, Bj, F, N), after the exclusion of the points lying more 
than three standard deviations from the mean. As most of the objects in the field are 
non-variable we assume as a measure of the r.m.s. photometric noise on the flux 
differences 5f Ui . Using the entire plate set, consisting of 5, 15, 5, 6 plates in the U, Bj, F, N 
bands respectively, it is also possible to derive a magnitude dependent photometric error 
for each plate (instead of considering the error on the magnitude difference in plate pairs). 
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As above we sorted the 2540 objects according to their Bj magnitude (on some reference 
plate), then for each object we considered the light curve, its time average and its mean 
square dispersion, in each band. Then we took 100 neighboring objects in magnitude and 
computed the r.m.s. dispersions around the time average. Neglecting the variance of the 
time average, this provides an estimate of the r.m.s. photometric error of individual plates. 
The errors estimated from plate pairs are consistent with a quadratic combination of the 
errors of the relevant individual plates. These errors are shown for each point of Figure 3. 

The variance of 5a is a linear combination of the variances of the signals in the bands 
involved, with appropriate coefficients depending on the frequencies of the bands. We can 
fit to the 5a,5logf Uj data a straight line taking into account the errors of individual points 
on both axes. The actual calculation is done using the "fitexy" subroutine from flPress &;' 



Tcukolsky 1992|) , which gives also the statistical uncertainties on a and b. The result, with 



one sigma errors, is a = (—8.49 ± 5.50)10 2 and b = (2.55 ± 0.75). 

Repeating the analysis with the exclusion of the N band reduces the effect of a 
deviation of the SED from a power low providing a less biased, but more noisy, estimate of 
the local spectral slope. The correlation is lower, but still positive and marginally significant 
(p = 0.27, P{> p) = 8.87- 10~ 2 ). The straight line fit to the data, a = (-1.67 ± 3.07)10~ 3 
and b = 2.29 ± 0.37, is quite consistent with the previous result. 



Di Clement e et al. 1996 found that the average increase of variability with 
redshift is A5i/Alog(l + z) ~ 0.25 — 0.30. This quantity is related to the slope b by: 
A5i/Alog(l +z) ~ ASi/Alogz/ ~ (5a 2 ) 1 / 2 ~ b ({5logf u ) 2 ) 1/2 , where (5a 2 ) 1 / 2 is the 



r.m.s. fluctuation of the spectral slope and Si is the variability indicator defined in |Dl 



Ulemente et al. 1996| and ((5 log f v ) 2 ) 1 / 2 — 0.1 is the r.m.s. flux density fluctuation. Thus 



the value of b found in the present analysis is consistent with the increase of variability with 
frequency found by |Di Clemente et al. 1996| and confirms the interpretation of the positive 
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v-z correlation suggested by GTV , by a direct observation of the r.m.s. slope variability 
We can compare our statistical results with the a — logf u relation found by |Edelson~e~t] 
|al. 19901 i n individual objects through the analysis of IUE observations. All of their 6 
objects have low redshift and are sampled around a wavelength A ~ 2000 A which is about 
the average wavelength of our sample (< z >= 1.365). For 5 of their objects they find a 
significant correlation of a and logf nu . For 4 of these objects the slope b of the relation 
a = b logf v + const is around 1.7, while 3C 273 shows a steeper {b ~ 2.8) relation. 

It would also be interesting to compare our findings with the variability of (B-R) colors 
measured by |Giveon et al. 199^| in a sample 42 of PG QSOs. Although the average redshift 
and the rest-frame wavelength region sampled, < z >~ 0.12, A ~ 3500 A differ from ours, 
qualitatively their results agree with ours, namely they also find a hardening of the spectrum 
during bright phases. A more detailed and quantitative comparison requires a re-analysis 
of |Giveon et al. 1999] data in a consistent way (Trevese & Vagnetti, in preparation; Trevese 



& Vagnetti 2000a,b). 



4. Discussion and Summary 

Let us assume that the SEDs of the objects of our sample are typical of QSOs, namely 
they are dominated by the big blue bump in the spectral region around A = 2000 A sec 
e.g. |Bregman 1990| ). We want to check the hypothesis that both the SED slope and 
brightness changes are caused by temperature changes of an emitting black body. For this 
purpose we use U, Bj, F data only and we assign the slope a to the intermediate frequency 
logu = \tog{yuVjVp). Figure 5 shows a versus logz/ for each object. In the same figure the 
curves a(x) = 3 — xe x /(e x — 1), x = hv/kT, of black bodies are also reported, for various 
values of the temperature T, h and k being the Planck and the Boltzmann constants. 
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EDITOR: PLACE FIGURE | HERE. 

The two points marked with stars represent the uncorrected SED of the two highest 
redshift objects, which are affected by the intergalactic Ly-a absorption and appear 
particularly steep. Since SA 57 is close to the Galactic pole, interstellar reddening has 
a negligible effect on the location of points in Figure 5. The distribution of a is skewed 
towards negative values, with most objects around -0.5. To check the consistency of 
slope and brightness variation in a more stringent way, we exclude from the sample the 
most deviant points, with a < —2a a . Then we obtain < a >= —0.48 ± 0.68. The 
corresponding average temperature is T m 25000i^. We cannot compare 5a/5logfp with 
a of individual objects because the noise is too large, especially on the former quantity. 
However we can compare the slope b of the straight line fitting the points in Figure 
4, representing the average increase of the spectral slope a for increasing luminosity, 
with the corresponding relation expected for black body spectra of varying temperature. 
Brightness and slope variations of the SED of a black body of fixed surface are related by 
{da/dT)/(dlogB v /dT) = (lnl0)[l - x/{e x - 1)] = F BB (x). In Figure 6, F BB (x) and a{x) 
are plotted as functions of x. 

EDITOR: PLACE FIGURE | HERE. 

The value of < a > and the r.m.s. spread a a of the sample define an interval of x in this 
plot. In the same figure is also reported the value of b and (b — er&) as deduced from the 
statistical analysis of section 2, which define a lower limit on x. 

As seen from Figure 6, the above values define a non-empty region of consistency 
between F BB (x) and a(x) for the sample. A black body with x ~ 3 satisfies simultaneously 
the two constraints. Thus we can say that the 5a and 5logf Uj , observed at an average 
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rest-frame A ~ 2000 A are consistent with temperature variations of a typical black body 
of temperature T « 25000K. 

A single black-body spectrum is an oversimplified representation of the SED which, in 
the case of accretion disks models, ignores temperature gradients and the presence of other 
components necessary to explain the emission outside the frequency range considered. In 
any case we can say that: 



i) The increase of the amplitude of variability with the rest-frame frequency (|Giallongo 



Trcvese fc Vagnetti 1991] , |Cristiani et al. 1996| , Pi Clemente et al. 1996|) is due to 



the hardening in the bright phases (and vice versa) of the spectrum, which to a first 
approximation maintains its local power-law shape. 

ii) The r.m.s. slope variations are quantitatively consistent with the increase of 
variability for increasing frequency found by |Di Clemente et al. 1996| . This confirms, by a 



direct observation of spectral slope changes, that correlation between the variability and 
redshift is accounted for by intrinsic spectral changes and does not require cosmological 
evolutionary effects. 

iii) Our results show that the average changes of the SED slope and the relevant flux 
variations in the optical- UV region are at least not inconsistent with temperature variations 
of a single emitting black body of T « 2.5 • 10 4 i^. 

iv) The average relation between slope and brightness variations provides a new 
constraint on models of the emission mechanism. 



For instance, an optically thin plasma model ( [Barvainis 1993| ) requires a temperature 
of the order of 10 6 i^ to explain with a single emission component the SED from optical-UV 
to X-ray frequencies (see [Fiore et al. 1995| ). The bremsstrahlung luminosity L u oc T~ x l 2 e x 
( |Rybicki fc Lightman 197^ ), can provide an optical-UV spectral slope a consistent with the 
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observations, but in this case the slope and brightness changes produced by temperature 
changes are related by Fff(x) = da/dlogf u — x/(x — 1/2), i.e. Fff ~ — 2x, for a; < 1. 
Thus, for T « 10 6 , Fff ~ —0.1 is not only too small in absolute value, but it is even 
negative, thus completely inconsistent with the result of our analysis. Notice that also 
the temperature changes of a simple black body of T « 10 5 K would be inconsistent with 
our findings since, for A ~ 2000 A, this temperature implies Fbb — 0.27. In fact, at high 
temperature any thermal emission approaches the limit where the long-wavelength part 
has a fixed spectral slope, independent of temperature changes ([Paltani and Courvoisici 



1994Q . Any model trying to reproduce the observed SED from the IR to X-rays could be 
tested against the observed a — logf u relation. In particular, it is possible to investigate 
whether the variation of one, or some, of the parameters involved (like the temperature) 
are capable of reproducing the observed relation. Perhaps additional phenomena, e.g. hot 
spots or flares ( [Krishan fc Wiita 1994] , |Kawaguchi et al. 1998|) over an otherwise stationary 



accretion disk, are needed to explain the SED variability. The main limit of the present 
analysis are the small amplitude of variability, due to the short time interval (1 yr) between 
the observations, and the lack of a — logf u relation for individual objects, which could be 
different depending, e.g., on their intrinsic luminosity. 

Despite these limits, our results clearly show that a few additional multi-band 
observations, properly distributed in time, could provide strong constraints on the physics 
of emission and variability, especially if associated with simultaneous observations in the IR 
and X-rays. 



This work was partly supported by the Italian Ministry for University and Research 
(MURST) under grant Cofin98-02-32 
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Fig. 1. — Redshift of all the objects of the sample vs. the apparent magnitude (m 4 ), defined 
as the Bj magnitude evaluated in a fixed aperture of radius of I'M. The stars represent 
the two highest redshift objects. The open circles represent the five objects with extended 
images from BTK. Lines of constant luminosity Mj = —23 and Mj = —21 are drawn 
assuming H Q = 50 km s _1 Mpc _1 , q Q = 0.5, anda = —1. 

Fig. 2. — U,Bj,F flux as a function of the rest-frame frequency, normalized at 2000 A, for 
the 40 objects of the sample. The photometric data of the 5 objects from the list of BTK98 
are represented by dots. 

Fig. 3. — The spectral energy distribution in arbitrary scale. Filled circles: first epoch; open 
circles: second epoch. The relevant power-law fits are also shown. The scale is the same for 
each pair of data sets representing the same objects at two different epochs. 

Fig. 4. — 5a versus 5logf v , with /„ measured in the Bj band, on independent plates (see 
text). The dashed lines represent the linear regressions, from which the correlation coefficient 
p is computed. The solid line is the linear fit which takes into the errors in both coordinates. 

Fig. 5. — The slope a as a function of x = hu/kT, for the 40 objects of the sample, 
as computed from U,Bj,F bands, with logv = \Ioq{vuVbVf)- The stars represent the 
two highest redshift objects whose colors are affected by Ly a absorption. The open circles 
represent the five objects with extended images from BTK. The solid lines represent a(x) 
for black bodies of different temperatures T, indicated in units of 10 4 K. The horizontal lines 
represent < a > ±cr Q , after the exclusion of the 3 points at more than 2-sigma from the 
mean. 

Fig. 6. — ct(x) and F BB (x) for a black body, as a function of x = hu/kT. The continuous 
horizontal lines represent respectively < a > and the slope b of the straight line fitting the 
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data in Figure 3. The dashed and dotted lines represent the one-sigma interval and the 
relevant bound on x, deriving from the comparison of the data with the black body curves 
Fbb(x) and a(x) respectively. 
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NOTE.— Ng94 is the serial number of point-like objects in T89, N2540 refers to the catalog described in BTK, N_, er is the serial number of the matching objects in the catalog of 
Koo 1996 as reported by Munn et al. 1997 and KKC is the serial number in Koo, Kron, & Cudworth 1986. The magnitude m.4 in column 7 is computed in a circular aperture 
of l".l on the Bj plate MPF1053 as in T89 and. The redshifts z are taken from T89 and T94. Columns from 9 to 17 contain the magnitudes in the same aperture of l".l of 
1114 obtained from the indicated plates. 
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